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ABSTRACT 


"V 

The  problem  of  detecting  high  altitude  rocket  plumes  with  satellite-borne 
submillimeter-wave  radiometers  is  examined  from  both  theoretical  and  experi¬ 
mental  points  of  view.  To  estimate  the  sizes  of  plume  signatures  contrasted 
against  a  250-K  earth  background  or  in  self-emission  against  the  cold  sky, 


a  computer  program  has  been  developed  to  predict  plume  brightness  temperatures 
and  optical  depths  of  rotational  lines  of  plume  molecular  constituents  (e.g., 
H^O)  as  a  function  of  distance  from  the  nozzle.  The  methods  employed  in  the 
computations  are  described  in  general  terms,  and  examples  are  presented  to  indi¬ 
cate  that  detectable  Hj,0  signatures  extending  to  several  thousand  nozzle 
diameters  should  exist  at  plume  altitudes  above  250  km.  Details  of  a  laboratory 
water  vapor  jet  operating  in  a  high  vacuum  environment  designed  to  simulate  a 
rocket  plume  at  high  altitude  are  outlined.  Sensitivity  considerations  relevant 
to  the  plume  detection  problem  are  discussed  for  Fourier  transform  spectrometers 
and  heterodyne  receivers  designed  to  study  spectral  lines  in  the  submillimeter- 
wave  region.  Finally,  a  technique  for  converting  the  heterodyne  radiometer  into 
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I.  INTRODUCTION 

The  possibility  of  observing  high  altitude  (i.e.,  above  150  km)  plumes 
of  third  stage,  mid-course  vernier,  or  bus  rocket  motors  from  satellites  by 
means  of  submillimeter-wave  heterodyne  detectors  was  discussed  in  an  earlier 
report  (l ).  Because  atmospheric  opacity  is  extremely  high  in  the  range 
of  frequencies  from  300  to  3000  GHz,  satellite  detection  of  rocket  launches 
would  not  be  disrupted  by  terrestrial  radiation,  either  natural  or  intended. 
Other  advantages  include  the  possibility  of  utilizing  agile  beam  electronic 
scanning  techniques,  the  likelihood  of  determining  altitude  and  velocity 
parameters  from  Doppler  shifts  of  spectral  lines,  and  the  benefits  associated 
with  the  expected  low-temperature  signature  of  the  plume. 

Through  the  approximately  adiabatic  expansion  which  occurs  as  the  rocket 
effluent  expands  from  the  nozzle,  rapid  cooling  creates  a  large  region  of  low 
temperature.  Surrounding  this  cold  core  is  an  envelope  of  higher  temperature, 
the  air-gas  mixing  layer,  created  by  collisions  between  plume  gas  and  atmos¬ 
pheric  molecules.  This  mixing  layer  region  is  believed  to  be  relatively  trans¬ 
parent  at  the  submillimeter  and  infrared  wavelengths  for  the  high-altitude  case. 
Since  some  of  the  low-energy  rotational  states  of  the  H20  molecule  ( a  prominent 
plume  constituent)  are  expected  to  remain  in  approximate  thermal  equilibrium 
with  the  adiabatic  expansion  of  the  core  region,  selected  spectral  lines  in 
this  frequency  range  should  be  optically  thick  and  colder  than  the  warm  earth 
background . 

Under  such  conditions,  radiometry  is  not  only  possible  within  the  width 
of  the  spectral  line,  but  offers  an  additional  advantage  of  not  being  dependent 
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on  the  existence  of  hot  portions  of  the  plume,  such  as  the  mixing  layer.  Any 
attempts  to  reduce  the  intensity  of  the  plume  infrared  radiation  by  modifying 
fuel  components,  for  example,  would  not  be  expected  to  decrease  the  size  of 
the  cold  signature  at  submillimeter  wavelengths. 

In  this  report,  the  problem  of  observing  cold  rocket  plumes  from  a 
satellite  is  examined  from  two  aspects:  1)  the  temperature,  absorptivity 
(discussed  in  terms  of  the  optical  depth  of  the  spectral  line),  and  size  of 
the  plume  signature,  and  2)  the  instrumentation  required  to  detect  a  plume 
both  in  laboratory  simulations  and  in  the  systems  application  described  pre¬ 
viously  (  1). 
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II.  THEORETICAL  STUDIES 

A.  Plume  Core  Versus  Mixing  Layer  Phenomenology 

To  describe  the  high-altitude  plume  signature  theoretically,  different 
sets  of  calculations  are  required  for  typical  streamlines  that  lie  along  the 
center  of  the  plume  core,  on  the  one  hand,  or  at  the  inner  edge  of  the  mixing 
layer  on  the  other  hand.  The  plume  core  is  that  part  of  the  rocket  exhaust 

gas  that  flows  downstream  and  is  initially  unaffected  by  collisions  with  the 

ambient  air.  The  mixing  layer  is  the  region  surrounding  the  plume  core  where 
ambient  air  and  plume  gas  mix  equally  and  interact  by  momentum- changing  colli¬ 
sions.  Eventually,  the  plume  core  is  dissipated  by  these  collisions,  and  the 
exhaust  gas  diffuses  into  the  atmosphere  a  few  mean-free  paths  downstream  (2). 

The  core  is  expected  to  be  radiationally- thick  in  the  transitions  between 
the  low-lying  rotational  energy  levels  considered  here.  The  mixing  layer  should 
be  radiationally- thin  in  these  same  transitions.  The  hot  mixing  layer  is 
caused  by  collisions  between  the  plume  gas  and  the  ambient  air  that  might  be 
moving  past  at  a  relative  velocity  of  about  5  km-s_1.  The  signal  above  the 
earth  background  originating  from  the  hot  mixing  layer  due  to  self -emission  is 
not  as  strong  as  the  signal  below  the  earth  background  from  cold  plume  core 
absorption  of  earth  background  radiation,  in  the  low-lying  rotational  transi¬ 
tions.  Differing  amounts  of  Doppler-shifting  in  mixing  layer  and  core  further 
distinguish  the  two  regions  in  a  frequency-resolved  spectrum  of  the  plume, 
whether  or  not  the  earth  is  viewed  as  background.  Plume  core  self-emission 
at  antenna  temperatures  ranging  between  about  50  and  100K,  against  the  cold 
(2.7  K)  sky,  is  also  expected. 
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The  basic  structure  of  a  rocket  plume  at  high  altitudes  is  shown  in 
Fig.  1,  where  the  inner  core  and  outer  mixing  layers  appear  as  distinct  regions. 
The  relative  velocity  of  the  two  fluids  relaxes  to  small  values,  com¬ 
parable  to  the  local  thermal  speed,  in  a  characteristic  mixing  layer  thickness 
(in  the  direction  of  given  by 


l 

m 


( vyv  )!l! lii 

^  7  12;  Pl+p2  a 


(1) 


where  m^,  and  m,,,  p2  are  the  mean  molecular  mass  (g/particle)  and  mass  den- 

sity  (g/cnr  )  of  the  ambient  air  and  plume  core,  respectively,  on  opposite  sides 

of  the  mixing  layer.  The  factor  a  is  the  effective  momentum- changing  collision 
2 

cross-section  (cm  ).  The  relative  velocity 


(2a) 


and  the  mean  fluid  velocity 


V  =  (p^  +  P2V2)/(P1  +  P2),  (2b) 

where  and  V2  are  the  fluid  velocities  on  opposite  sides  of  the  mixing 
layer  at  a  point  downstream. 

The  temperature  of  the  mixing  layer  is  calculated  from  the  specific 
enthalpy  h  of  the  mixture,  whose  stagnation  value  is  the  mass  average  of  the 
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EARTH  BACKGROUND 


Fig.l.  Core  and  mixing  layer  regions  of  high  altitude  rocket  plume 
radiometric  signature. 


stagnation  values  of  specific  enthalpies,  h^  and  h2,  of  the  two  fluids  on 
opposite  sides  of  the  mixing  layer.  Therefore,  making  the  approximation  that 
the  mean  fluid  velocity  V  is  nearly  constant  across  the  width  of  the  mixing 
layer  at  a  given  point  downstream,  and  assuming  that  the  flow  is  stationary 
in  the  coordinate  system  fixed  in  the  vehicle,  we  obtain 

h  *  |  v2  ♦  .  y  (hj  *  |vf)  *  (l-y^thj  ♦  |v§)  (3a) 

P 

where  V  2  is  the  local  velocity  difference  within  the  mixing  layer,  where  the 
local  fraction  of  atmospheric  species  has  the  value  y^  =  p^/p,  and  p  =  p^+p2 
is  the  total  mass  density.  Note  that 


Y1  \  +  (l-yi)  |  V2  =  |  V2  +  y  (l-y  )  j  V^, 


(3b) 


where  V^2  =  |V^-  V2| .  Subscript  1  refers  to  the  atmosphere  and  2  refers  to 
the  plume. 

The  mixture  temperature  T  is  a  maximum  T^  within  the  mixing  layer  when 

% 

the  local  relative  velocity  V^2  is  zero  and  when  the  mass  fraction  of  en¬ 
trained  air  y^  is  given  by 


y 
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1 
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(4) 
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The  approximation  of  constant  specific  heats  may  be  used  downstream,  so  that 
hcsiCpT,  and  ^2~Cp2^2*  Therefore,  the  temperature  maximum  is  approxi¬ 

mately 


T  ~ 
m 


c  ,T  + 
Pi  1 


c  0T_ 

P2  2 


i  V2 

4  12 


2  +  C 

pi  p2 


+  c 


'pi  p2 


(5) 


For  the  following  parameter  values, 


V12  =  5  km/s 

=  k/n^  =  0.83  *  108  erg/K-g 
c  ^  =  0.14  (64?  02,  36?  N2  by  volume) 

cp2  =  0.19  (30%  H2,  30%  CO  or  N2,  10?  C02,  30%  H20  by  volume) 

Tp  =  1000  K 
T2  =  100  K 

for  which 


Cp  =  yi  Cpl  +  (l"yl}  °P2  =  °-165  *R 

we  obtain  T  =  2750  K,  of  which  only  480  K  arises  from  the  firkt  term  of 
m 

Eq.  (5). 

The  ratio  of  water  vapor  optical  depths  for  the  mixing  layer  and  the 
core  for  a  given  radiative  transition  is  approximately  given  by 
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(6) 


_  O 

where  n^  and  n^  are  the  water  vapor  concentrations  (cm  )  in  the  mixing  layer 

and  the  core,  and  M  is  the  local  Mach  number  in  the  core  at  this  distance  r 
c 

downstream.  The  temperatures  with  subscript  x  denote  the  excitation  tempera¬ 
tures  for  the  two  levels  that  yield  the  fractional  population  in  the  lower  of 
the  two  levels  (separated  by  energy  E);  and  with  subscript  t,  the  translation 
temperatures  that  determine  the  thermal  speeds. 

To  obtain  a  rough  estimate  of  the  ratio  of  optical  depths,  we  set  all 
three  core  temperatures  to  100  K  and  all  three  mixing  layer  temperatures  to 
1000  K.  The  Boltzmann  factor  in  Eq.  (6)  is  approximated  by  the  value  2  for 

the  557-  and  752-GHz  lines,  and  the  Mach  number  M  ^  4.  For  a  nozzle  diameter 

c 

d,  the  core  H20  density  (in  the  absence  of  i^O-air  collisions)  is  approximately 


,-.17,..  x2  -3 
n  ~  10  ( d/r )  cm 

o 


(7a) 


The  mixing-layer  H^O  density  is  approximately 

-t/Z  n  r 

n  =  n  -  n  =  n  (l-e  C)  —  for  t/Z  <  1 ,  (7b) 

m  o  c  o  Z  c 
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-r/Jl 

where  =  nQe  c  is  the  H^O  core  density  in  the  presence  of  H^O-air 

collisions.  Beyond  r  —  &Q>  nc  decreases  exponentially  so  that  the  core  is 

no  longer  important.  The  mean-free-path  for  H^O-air  collisions, 

-15  2 

Jl  =  m^/(a'p^)~l  km  at  250-km  altitude  with  a'^slO  cm  as  the  H^O-air 
collision  cross-section. 

We  then  obtain  from  Eq.  (6)  the  ratio  of  optical  depths 


T  q 

—  ^  4  *  10^  for  r/SL  <  1. 

T  C 

C 


(8) 


Assuming  that  we  are  calculating  the  ratio  of  antenna  temperatures  for  the  mix¬ 
ing  layer  and  core  sufficiently  downstream  that  the  core  is  just  thin,  then 


AT  (T  -T)  t 

I2TT  =  rrbrT  f  ~  for  r/1c  <  (9) 

c  c  c 

where  T  is  the  earth  background  temperature.  Thus,  we  conclude  that,  though 
the  mixing  layer  is  hot,  it  is  too  thin  to  affect  the  signal  due  to  the 
core,  even  if  there  were  no  Doppler  shift  between  the  two  regions. 

B.  Plume  Signature  Calculations 

A  computer  program  (Plume  1)  has  been  developed  for  the  MIT  Lincoln 
Lab  IBM/370  machine  to  predict  the  intensities  of  the  rotational  lines  of 
plume  molecules  such  as  water  vapor  or  carbon  monoxide  at  various  distances 
from  the  rocket  nozzle. 
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For  the  present,  simplifying  assumptions  limit  the  results  of  the  calcu¬ 
lations  to  the  high-altitude  case  (above  150  km),  where  no  shock  phenomena  are 
expected  to  occur  downstream  because  the  interaction  between  plume  gas  and  ambient 
air  is  too  diffuse.  The  air  mean-free-path  for  momentum-changing  collisions  is 
at  least  100  meters  long  at  these  altitudes. 

For  water  vapor,  the  radiative  transport  is  calculated  for  all  electric- 
dipole  transitions,  however  weak,  between  pairs  of  rotational  energy  levels 

up  to  angular  momentum  quantum  number  J  =10.  A  value  of  J  =10  has  been  found 

m  m 

to  be  just  large  enough  to  represent  the  bulk  of  the  distribution  of  populations 

among  the  rotational  levels  for  translational  temperatures  less  than  1000  K  that 

are  encountered  downstream.  There  are  605  electric-dipole  radiative  transitions 

for  J  =10.  In  general,  the  number  of  transitions  N  for  an  asymmetric  rotor 
m  s 

2  2 

(e.g.,  H  O)  is  given  by  N  =(l/2)J  N  ,  where  N=(J  +1)  is  the  number  of  energy 
c  s  m  m 

levels.  This  cubic  and  quadratic  dependence  on  contrasts  with  the  much  sim¬ 
pler  case  of  the  linear  (diatomic)  molecule,  carbon  monoxide,  for  which 

N=J  +1  and  N  =J  ,  where  both  numbers  rise  only  linearly  with  J  . 
m  s  m  J  J  m 

The  line  intensities  are  calculated  at  present  under  the  simplifying 
assumption  of  the  "local  approximation."  This  approximation  is  most  useful 
where  supersonic  flow  is  present  because  a  photon  is  emitted  at  a  point  in  the 
flow  with  a  resonant  frequency  that  lies  within  a  narrow  range  about  the 
Doppler  shift  derived  from  the  component  of  the  mean  flow  velocity  along  the 
photon  propagation  direction. 

After  the  photon  propagates  over  a  short  distance,  its  frequency  will  no 
longer  be  resonant  in  this  transition  with  the  molecules  at  this  new  location 


because  they  are  moving  with  a  different  velocity  component  that  is  outside 
the  narrow  range  required  for  resonance.  The  narrow  frequency  width  is  the 
result  of  small  thermal  (or  micro-turbulent )  velocities. 

The  flow  is  approximatly  radially-directed  from  a  point-like  source  at 
the  center  of  the  nozzle  orifice.  Generally,  the  vehicle  velocity  is  not 
aligned  axially,  i.e.,  not  collinearly  with  the  central  exhaust  streamline. 

The  streamlines  are  then  pushed  sideways  by  the  broadside  momentum  transfer 
between  air  and  plume  streams.  Figure  2  shows  the  plume  core  and  mixing  layer 
regions  of  the  rocket  exhaust  for  such  a  misaligned  case. 

The  following  sections  describe,  in  general  terms,  how  the  thermodynamic, 
radiative  and  collisional  properties  were  incorporated  into  the  program. 

1.  Thermodynamic  Properties 

To  demonstrate  the  capabilities  of  the  program,  the  altitude  range 

of  200  to  250  km  was  chosen  for  most  of  the  calculations.  This  corresponds  to 

9  -3 

an  average  air  density  of  about  4  x  10  particles-cm  ,  and  a  pressure  of 
4  *  10~10  atms. 

The  density  and  temperature  of  plume  gas  emerging  from  the  nozzle  were 
17  -3 

selected  as  10  cm  and  1000  K,  respectively,  and  are  representative  of  most 

2  -1 

engines.  The  gas  density  was  taken  to  vary  downstream  as  (l+R  )  ,  where  R 

is  the  ratio  of  radial  distance  downstream  to  the  characteristic  distance  over 
which  the  density  falls  to  half  the  exit  value.  On  the  central  axis  this  dis¬ 
tance  is  slightly  larger  than  the  nozzle  diameter  in  the  Brook  model(3).  De¬ 
tails  of  the  flow  and  the  transverse  dependence  of  properties  through  the 
mixing  layer  will  be  described  in  another  report. 


Fig. 2.  Model  of  plume  core  and  mixing  layer  radiometric  signals  for 
high  altitude  plumes  as  a  function  of  frequency  for  arbitrary  line 
of  sight. 


The  translational  temperature  of  the  plume  gas  is  obtained  from  the 
energy  equation  that  sets  the  rate  of  energy  transfer  from  the  air  and  from 
the  internal  degrees-of-freedom  of  the  plume  gas  to  the  translational 
degrees-of-freedom  of  the  plume  gas  equal  to  the  sum  of  the  rate  of  increase 
of  the  translational  thermal  energy  and  the  rate  of  work  done  by  volume  ex¬ 
pansion.  To  compute  the  translational  temperature,  the  rate  coefficients  for 
transfer  of  translational  energy  from  atmosphere  to  plume  gases  should  range 
between  1  and  10  s-1  for  the  200  to  250-km  altitude  range. 

2.  Radiative  Properties 

The  line  intensity  and  line  shape  of  any  transition  is  calculated 
as  an  excess  radiation  over  the  background  continuum  radiation  corresponding 
to  the  temperature  of  2.7  K  of  the  cold  sky  or  250  K  of  the  earth's  atmosphere. 

The  intensity  AI^v  of  the  excess  radiation  over  the  background  is  an 
integral,  with  respect  to  optical  depth  of  the  product  of  exp(-x^v)  and 

the  difference  between  the  source  function  S  and  the  background  intensity  1^. 
The  optical  depth  x^  is  the  integral  with  respect  to  distance  along  a  ray 
path  (with  direction  cosine  y)  of  the  absorption  coefficient  at  frequency  v. 

AI  =  /  dx  e  yv  ( S-I.  ) .  (10) 

yv  1  yv  b 

In  the  local  approximation,  the  optical  depth  is  a  quantity  associated 
with  each  point  in  the  flow  and  associated  with  the  Doppler-shifted  resonance 
frequency  obtained  from  the  component  of  the  flow  velocity  along  the  ray 
direction. 
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The  source  function  is  the  intensity  given  by  the  Planck  function 


evaluated  for  the  excitation  temperature  T  ,  which  is  the  effective  Boltzmann 
temperature  between  the  pair  of  levels  whose  transition  is  being  considered. 

The  excitation  temperature  is  inversely  proportional  to  the  logarithm  of  the 
ratio  of  the  two  sublevel  populations.  Thus,  the  solution  of  the  rate  equations 
for  the  populations  provides  both  T^  and  x  at  each  point  in  the  flow  for 
each  transition. 

When  the  local  approximation  is  applied  to  the  integral  expression  for 
the  excess  intensity,  each  point  in  the  flow  and  each  direction  through  that 
point  has  its  assigned  intensity  obtained  from  the  formula 


AI 


Viv 


2hv 

7* 


tivAT-x  , 
e  -1 


1 

hvATb 

e  -1 


( 1-e 


(11) 


where  Tfe  denotes  the  temperature  of  the  background.  In  terms  of  antenna  temper¬ 
ature  T^,  the  effective  temperature  difference  between  plume  and  background, 

Eq.  (ll)  becomes 


T.  =  AI  X2/2k  .  (12) 

A  pv  ' 

The  local  approximation  to  the  integral  in  Eq.  (10)  need  not  be  used  if 
the  line  shape  is  required  to  reflect  accurately  the  Doppler  effects  of  thermal 
and  turbulent  motions,  as  well  as  those  of  the  mean  flow.  This  approximation 
might  still  be  retained  by  way  of  the  photon  escape  probabilities  to  obtain 
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( from  the  approximate  rate  equations )  the  population  densities  that  appear 
in  the  excitation  temperatures  and  the  optical  depths  which  are  to  be  used 
in  the  integral  for  the  intensity. 

3.  Collisional  Excitation  Properties 

The  transition  rate,  due  to  thermal  collisions,  between  any  pair 
of  energy  levels  separated  by  energy  AE  is  given  by  the  phenomenological  formula 
for  the  de-excitation  rate  coefficient  (s-^  per  molecule) 

Ci  =  Cj  (1  +  a^x )  exp(-b.x),  (13) 

where  i  labels  the  first  of  the  colliding  pair  (e.g.,  H20  -  H^,  0  -  H20,  etc.) 
and  x  =  AE/M\.  The  temperature  T^  is  the  translational  temperature  of  the 
species  colliding  with  ^0,  if  H20  is  the  radiating  species  being  considered 
in  the  rate  equation.  If  atomic  oxygen  of  the  ambient  atmosphere  is  causing 
the  excitation,  then  T^  =  Tq  =  1000  K.  In  this  case,  the  parameters 
Sq  =  1  and  b^  =  2  were  chosen.  For  -  H^O  collisions,  T^  =  T,  the  plume 

translational  temperature  at  the  point  in  question,  and  a^  q  =  b^  ^  =  1  were 
chosen.  These  values  of  the  a  and  b  parameters  are  selected  based  on  previous 
experience  with  theoretical  and  experimental  values  known  for  the  collision 
pairs  H2  -  CO  and  HF  -  HF.  It  should  be  noted  that  if  a^  is  greater  than  b^, 
the  rate  shows  a  maximum  as  a  function  of  AE,  which  is  a  quantity  that  in¬ 
creases  for  two  levels  directly  as  the  difference  between  their  values  of  J, 
the  total  angular  momentum  quantum  number.  For  the  conditions  under  con- 
sideration,  the  value  of  C!  for  plume  molecules  is  10  cm  -s  ,  and  for  both 
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entrained  air  at  temperature  T  and  ambient  air  at  Tq,  the  value  of  Cj  is 
-11  3  -1 

3  x  10  cm  -s  .  By  comparison  with  this  rate  constant  for  rotational 
excitation,  the  rate  constant  for  heating  of  the  plume  gas  translational 
temperature  by  the  air  collisions  at  temperature  Tq  is  taken  to  be  ten  times 
larger. 

High-speed  collisions  also  occur  because  of  the  relative  flow  velocity 
that  exists  between  the  flow  of  plume  gas  and  the  ambient  air.  The  inelastic 
collision  rate  coefficient  for  deexcitation  of  rotational  levels  in  the  plume 
gas,  owing  to  this  directed  flux,  is  of  the  form 

C  =  C'/F  s"1  (14) 

and  the  excitation  rate  coefficient  is 

C  =  C'(l-xom)/F  (X0„<1)  (15) 

c  =  o 

where  x  =  AE/E  and  E  is  the  energy  available  in  the  center-of-mass  system 
of  a  collision,  while  F  is  a  factor  that  makes  the  collisional  excitation  of 
H-,0  levels  with  lower  values  of  quantum  number  K  more  probable.  This  trend 
was  noted  from  quanta!  calculation.;  made  earlier  for  inelastic  collisions  of 
0  and  HpO  (internal  Lincoln  Laboratory  memorandum).  For  these  nigh-speed 
collisions,  C'  is  chosen  to  be  10  ^  cm^  -s  ^  (for  200-250  km  altitudes  and 
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a  relative  velocity  of  5  km-s-"'").  For  the  plume  core  (near  the  inner  edge  of 
the  mixing  layer),  only  10%  of  the  directed  flux  is  assumed  to  be  participating 
and  the  effective  C*  becomes  10-^  em^-s-\ 

4.  Photon  Escape  Probabilities 

The  factor  exp( -x  )  is  the  probability  that  a  photon  of  frequency  v 
will  travel  a  distance  corresponding  to  the  optical  depth  in  the  direction 
specified  by  the  direction-cosine  p.  When  averaged  over  solid  angle  at  a  point 
and  over  frequency,  as  weighted  with  the  local  frequency  distribution  function 
f(v-v  )  for  thermal  and  microturbulent  velocities,  this  factor  becomes  the 
photon  escape  probability  B  that  multiplies  the  Einstein  A-  or  B-coefficient 
in  the  appropriate  rate  equation.  The  probability  modifies  the  radiation 
processes  of  spontaneous  and  stimulated  emission  and  absorption  for  each 
transition.  The  frequency  distribution  function  f(v-VQ)  is  taken  to  be  a 
Gaussian  function  centered  on  the  resonant  frequency  Vq  that  is  Doppler-shifted 
by  the  parallel  component  of  the  mean  flow  velocity  and  Doppler-broadened  by 
the  root-mean-square  random  speed  due  to  motion.  The  distribution  function  is 
normalized  to  have  unity  for  its  integral  over  frequency. 

In  the  local  approximation  and  the  case  of  radial  source  flow  with 
velocity  proportional  to  the  radius  from  the  point  source,  the  photon  escape 
probability  is  given  by 
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B  =  (l-e"T)/T, 


where 


T  =  otr/v 


Ag  n.  n  0  , 

—  &u  /  £  u\  -,3  -1 

a  =  {ST-  (- - — *  S 

8,7  gu 


is  the  absorption  coefficient  after  integration  over  frequency  and  after  multi¬ 
plying  by  wavelength  X  to  convert  to  Doppler-velocity  units  for  frequency- 
offset  from  line  center.  The  population  densities  and  degeneracies  in  the 
lower  and  upper  levels  of  the  transition  are  denoted  by  n^  and  n^,  and  g^  and 
gu  respectively. 

For  the  case  of  constant  velocity  in  point-source  flow,  the  photon 
escape  probability  becomes 


0^-e 


[K2^  "  "  K1^2^* 


This  dependence  of  3  on  modified  Bessel  functions  of  the  second  kind  can  be 
shown  to  behave  sufficiently  like  the  exponential-type  function  of  the  pre¬ 
vious  case  that  the  exponential  version  is  used  for  the  present,  even  though 
the  constant  velocity  approximation  is  more  realistic  over  the  portion  of  plume 
of  interest  (i.e.,  where  the  temperature  is  low  compared  to  the  stagnation 
temperature ) . 
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The  A-coefficient  is  determined  from  the  relative  line  strength  as  cal¬ 
culated  from  the  computer  program  TOP,  provided  by  Harvard  University  Chemistry 
Department,  for  asymmetric  rotors  with  centrifugal  distortion.  The  rotational 
energies,  however,  are  calculated  from  a  program  developed  by  Cook  et  al.(4) 
and  provided  by  Prof.  Frank  Delucia  of  Duke  University.  This  program  corrects 
for  centrifugal  distortion  more  accurately  so  that  the  line  frequencies  are 
accurate  to  better  than  1  MHz  up  to  J=10. 

5.  Computational  Method 

The  rate  equation  for  the  fractional  population  density  n.  in  the 
ith  energy  level  has  the  form 

3n. 

V5F-'  BJn(l’r)  '  Rout(1'r)’  U8> 


where  R  , ,  the  rate  of  departure  from  this  level,  is  given  by 
out 


Eout<1’r)  ■  Vi- 


(19) 


The  coefficient  is  the  sum  of  all  rate  coefficients,  radiative  and  colli sional, 
that  carry  population  from  the  ith  level  to  any  other  level.  The  flow  velocity 
is  v,  which  is  a  function  of  distance  through  the  Bernoulli  equation  for 
diabatic  conditions.  The  solution  to  the  above  equation,  treated  formally 
as  linear,  is 
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-Oi(ro,r)  r  -ai(r',r) 

n.(r)  =  n.(r)e  +  f  ■,  e  R.  (i,rf),  (20) 

1  l  o  ■*  v(r '  )  m 

r 

o 

where  a.(r',r)  =  /  — g^(r' )•  The  exponential  factors  are  the  probabilities 

that  population  sources  at  upstream  points  r'  of  a  streamline  are  felt  down¬ 
stream  at  the  point  r.  For  the  present,  straight  radial  streamlines  may  be 
imagined,  but  these  are  bent  by  interaction  with  the  ambient  air  stream  accord¬ 
ing  to  flow  calculations  to  be  described  in  another  report. 

The  source  term  R.  is  the  sum  of  all  sources  from  all  levels  denoted 
in 

by  the  index  j,  that  is, 


R  (i,r’)  =  E.n.(r')  W..(r'), 

J  J  J 

where  W..  =  A..(l+I..)  g.  .  +  C..  (if  E.  >  E. ) 

Ji  J1  1J  1J  J1  J  i 


(21) 


A.  .(g./g.)  I . .  3  . .  +  C  . . 

i  J  Ji  Jl  Ji 


(if  E.  <  E. ), 
J  i 


depending  upon  whether  E.  is  greater  or  smaller  than  E^,  the  energy  of  the 
ith  level.  An  energy  level  diagram  for  the  rotational  levels  of  H^O  is  shown 
in  Fig.  3. 

The  non-dimensional  intensity  I.  .  that  appears  above  is  the  incident  or 

background  intensity  associated  with  the  earth's  atmosphere  (or  the  nozzle 

2 

orifice)  divided  by  2hv/A  .  The  collision  coefficient  C. .  (and  C..)  is  the 
'  J  l.i  ,n 


20 


ENERGY  (cm- 


appropriate  sum  of  the  rate  coefficients  for  all  the  different  collision  part¬ 
ners  that  induce  transitions  between  levels  i  and  j  in  the  radiating  species, 
e.g.,  H20. 

The  equations  for  the  n^(r)  are  solved  by  iteration,  i.e.,  by  introducing 
the  previously-obtained  n^(r' )•  The  initial  guesses  are  those  of  thermal 
equilibrium  at  the  local  translational  temperature.  For  the  present,  all 
numerical  integrations  are  done  with  the  simple  trapezoidal  rule  on  a  grid  of 
points  in  the  radial  position  variable  R,  the  ratio  of  the  position  r  to  the 
characteristic  length  in  the  density  law  discussed  earlier.  The  grid  spacing 
is  made  wider  for  points  farther  downstream,  where  properties  are  found  to  vary 
more  slowly. 

6.  Plume  Signature 

For  the  work  reported  here,  the  case  of  the  streamlines  of  the 
plume  core  that  lie  near  the  inner  edge  of  the  mixing  layer  will  be  discussed. 
As  shown  in  Fig.  4,  the  optical  depth  for  the  557-GHz  transition  in  water 
vapor  may  be  substantial  at  distances  downstream  as  great  as  4000  nozzle 
diameters,  for  an  altitude  of  250  km,  and  15,000  nozzle  diameters  for  350  km. 
Within  this  distance,  the  absorption  of  the  background  earth  radiation  would 
be  significant  if  the  plume  were  viewed  from  above.  In  Fig.  5,  the  excitation 
temperatures  for  the  557-GHz  and  752-GHz  transitions  rise  farther  downstream 
(see  Section  2),  due  to  collisional  heating  by  the  ambient  air  stream,  sc  that 
both  the  optical  depths  and  differences  between  source  function  and  background 
intensity  are  reduced.  The  diabatic  temperature  is  the  translational  plume 
gas  temperature  taking  into  account  heating  by  air-interaction  and  cooling  by 
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R  (Nozzle  Diameters  Downstream) 

Fig. 4.  Optical  depths  at  250  km  for  557-  and  752-GHz  water  lines 
as  a  function  of  axial  distance  from  the  exhaust  nozzle. 


T  (K) 


volume  expansion.  The  adiabatic  temperature  indicates  this  cooling  effect 
when  the  air  collisions  are  absent.  The  quantity  v/d  is  the  ratio  of  local 
exhaust  velocity  in  the  core  (nearly  constant  with  distance  downstream)  to 
the  nozzle  diameter.  This  quantity  characterizes  a  typical  flow  rate  (s-''") 
to  be  compared  with  collisional  and  radiative  rates. 

Figure  6a  shows  the  resulting  behavior  of  the  line  intensities  along  a 
typical  streamline  (given  in  terms  of  antenna  temperature  T^)  for  the  emerging 
^0  radiation  at  the  557-  and  752-GHz  frequencies.  In  Fig.  6b,  the  antenna 
temperatures  of  the  557-GHz  line  are  compared  at  250-  and  350-tan  altitude, 
where  it  is  seen  how  the  plume  signature  can  increase  dramatically  as  the 
ambient  air  density  falls  off  at  even  higher  altitudes.  In  Fig.  7,  an  esti¬ 
mated  typical  557-GHz  plume  signature  is  shown  as  a  function  of  frequency 
offset  (from  rest  frequency)  against  the  absorption  profile  of  the  earth's 
atmosphere.  The  computer  program  is  being  extended  to  yield  more  accurate 
radiation  intensities  and  profiles  from  detailed  models  of  the  mixing  layers 
and  core. 

From  these  calculations,  it  may  be  concluded  that  rocket  motors  with 
thrusts  of  only  a  few  thousand  pounds  at  altitudes  of  250  tan  may  produce  cold 
plume  signatures  large  enough  to  provide  sufficient  beam  fill  for  observation 
by  submillimeter-wave  heterodyne  radiometers  in  synchronous  orbit.  At  higher 
altitudes  (e.g.,  above  350  km),  motors  with  thrusts  of  a  few  hundred  pounds 
may  be  detectable.  Even  with  present  far-infrared  technology,  observation  of 
high-altitude  plumes  may  be  accomplished  with  antenna  apertures  of  several 
meters  in  diameter  and  dwell  times  on  the  order  of  one  second. 
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Fig. 6a.  Plume  radiometric  signals  for  557-  and  752-GHz  at  250-km 
altitude  against  a  250-K  earth  background  as  a  function  of  axial 
distance  from  the  exhaust  nozzle. 
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Fig. 6b.  Plume  radiometric  signals  at  557  GHz  against  a  250-K  earth 
background  as  a  function  of  axial  distance  from  the  exhaust  nozzle, 
at  altitudes  of  250  and  350  km. 
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III.  EXPERIMENTAL  TECHNIQUES 


To  examine  experimentally  the  theoretical  predictions  of  plume  rotational 
temperatures  and  optical  depths  of  submillimeter-wave  absorption  lines,  two 
basic  approaches  are  available  for  use  with  a  laboratory-scale  simulated 
plume.  For  surveying  the  entire  spectral  region  of  interest,  a  Fourier 
transform  spectrometer  is  potentially  the  most  appropriate  instrument,  and 
for  studying  the  characteristics  of  selected  absorption  lines,  a  heterodyne 
receiver  can  provide  the  necessary  capabilities.  In  both  cases,  sensitivity 
requirements  for  plume  applications  reach  to  current  state-of-the-art  limi¬ 
tations  . 

A.  Simulated  Water  Vapor  Plume 

To  diagnose  the  absorption  (or  emission)  characteristics  of  an 
abiabatically  expanding  rocket  plume  with  high  water  vapor  content,  a  large 
high-vacuum  chamber  has  been  adapted  to  accommodate  the  operation  of  a 

laboratory-scale  water  vapor  jet.  Depicted  in  Fig.  8,  the  Lincoln  Laboratory 
space  environmental  test  chamber  is  approximately  cylindrical  with  a  diameter 
of  three  meters  and  a  length  of  ten  meters.  Since  the  inner  surfaces  are 
almost  entirely  liquid-nitrogen  cryopaneled,  the  chamber  is  idealxy  suited 
for  pumping  water  vapor  and  is  capable  of  reaching  a  background  pressure  of 

_7 

3  x  10  Torr  (simulating  an  altitude  of  250  km).  Since  the  vapor  pressure 
of  ice  is  virtually  non-existent  at  77  K,  the  jet  effluent  can  continue  to 
freeze  on  the  chamber  walls  without  an  appreciable  rise  in  pressure  for 
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several  hours  at  mass  flow  rates  on  the  order  of  1  g/sec.  In  this  manner,  the 
high  altitude  conditions  of  a  rocket  plume  core  may  be  simulated  with  good 
accuracy. 

Figure  9  is  a  sketch  of  the  transparent  gold  furnace  (5  )  used  to  gen¬ 
erate  the  water  vapor  jet.  Because  of  the  absence  of  out-gassing,  these 
furnaces  are  ideal  for  in-vacuum  applications,  where  only  an  infrared 
radiation  shield  is  required  as  insulation.  A  thin  gold  film  evaporated  onto 
the  internal  surface  of  the  outer  pyrex  tube  provides  excellent  reflective 
properties  in  the  infrared,  while  rendering  the  furnace  transparent  at  visible 
wavelengths.  The  inner  reaction  tube  is  made  of  quartz  and  is  divided  into 
two  stages,  with  a  baffle  as  a  water  dam.  In  the  first  stage,  called  the 
vaporization  section,  the  water  is  brought  to  boiling.  With  the  escaping  steam 
passing  over  to  the  steam  section,  recommended  stagnation  conditions  of  900  K 
temperature  and  10  Torr  pressure  (6  )  may  be  achieved  with  modest  amounts  of 
heating  power.  The  dimensions  of  the  simulated  plumes  may  be  tailored  accord¬ 
ing  to  nozzle  diameters,  and  where  necessary,  by  increasing  the  size  and  heat¬ 
ing  capacity  of  the  gold  furnace. 

B.  Fourier  Transform  Spectroscopy 

As  discussed  in  the  Appendix,  incoherent  detection  by  means  of  a  gratiiig 
spectrometer  is  hampered  by  throughput  limitations  where  high  resolution  is 
desired.  A  Fourier  Transform  Spectrometer  (FTS)  has  the  advantage  of  a  ihrough 
put  that  is  independent  of  resolution  because  the  element  width  Av  is  determine 
by  the  maximum  distance  of  mirror  travel  L  (i.e.,  Av  =  1/.1L).  Since  the  ins  +  ru 
ment  is  essentially  a  Michelson  interferometer  (see  Fig.  8)  which  continually 
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Fig. 9.  Laboratory-scale  water  vapor  jet  with  sonic 


collects  all  of  the  spectral  radiation  from  a  given  collimated  source  over  a 
large  predetection  bandwidth  ( determined  by  the  spectral  band  of  the  detector 
and  any  filters  or  windows  in  the  optical  path)  as  the  moving  mirror  scans 
through  a  distance  L,  it  is  ideally  suited  for  survey  applications.  This 
multiplexing  advantage  means  that  an  FTS  can  scan  with  equivalent  sensitivity 
the  entire  spectral  region  in  the  same  time  that  a  grating  instrument  would 
require  to  cover  a  single  resolution  element  Av. 

For  detector-noise  limited  operation  in  the  Rayleigh- Jeans  region  (where 
hv  «  kT  ,  IcT^ ) ,  Eq.  (A-4)  and  Eq.  (A-9)  for  the  ratio  of  signal  to  rms  noise 
may  be  reduced  to 
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(§)  =  n1 - 2 - - - 2 —  (22) 

FTS  c  ( NEP). 

me 

where  AT  =  (T  -T,_)(l-e_T)  and  t  is  the  total  time  of  observation  (or  inte- 
P  b 

gration).  Negative  values  of  AT  denote  absorption  and  positive  values  denote 
emission.  The  optical  depth  t  of  the  plume  appears  in  the  exponential  function. 
For  this  type  of  instrument,  the  collection  efficiency  is  determined  by  the 
beamsplitter  loss  and  any  other  losses  associated  with  windows,  filters  and 
lightguides;  the  observation  efficiency  is  fixed  by  the  chopping  cycle 
associated  with  synchronous  detection  and  is  normally  equal  to  1/2.  For  a 
step-and-integrate  FTS,  t  represents  the  product  of  the  sampling  time  (time 
for  each  step)  and  the  total  number  of  steps  in  one  mirror  scan.  For  the  con¬ 
tinuous-scan  version,  t  is  the  time  for  one  scan,  multiplied  by  the  number  of 
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scans.  Since  mirror  velocities  are  much  greater  in  continuous  scan  models, 
there  is  usually  time  to  carry  out  several  scans  in  one  experiment 

With  the  particular  instrument  in  use  at  Lincoln  Laboratory  for  the 
spectroscopy  measurements  of  an  optically  thick  (  t  >  1 )  plume  with 

=  50  K,  for  example,  the  values  of  the  experimental  parameters  are  as 
follows : 

=  3000  K  (Hg  arc  lamp) 

6v  =  10"5u  (Ref.  2) 

A£!  =  1.6  cm^-ster  (23) 

(NEP)^.^  =  5  x  10  13  W-Hz  (Ge  bolometer) 

n-jL  <  0.25 

n2  <  0.5 

Thus,  Eq.  (22)  may  be  simplified  to 

(§)  <  5  x  io“7  v3^,  (24) 

FTS 

with  v  in  GHz.  For  v  =  752  GHz  and  a  total  integration  time  t  =  1  s  (achievaDle 
with  a  rapid-scan  instrument) 


(£) 

^  N 


<  212 


FTS 


(25) 
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for  the  752-GHz  water  vapor  absorption  line  (2q2  -*  2^,  see  Fig.  3).  Similar 
results  apply  for  the  557-GHz  line  (1  -*■  I^q). 

The  above  estimate  of  the  spectrometer  sensitivity  suggests  that  absorp¬ 
tion  measurements  of  plume  lines  may  be  carried  out  with  little  difficulty. 
However,  this  computation  is  based  on  the  assumption  that  the  system  is 
limited  to  detector  noise  only  and  that  the  full  multiplexing  advantage  of  the 
FTS  can  be  realized.  The  extent  to  which  this  assumption  is  valid  depends  on 
the  intensity  of  other  noise  sources  which  may  dominate  the  detector  noise 
contribution.  For  a  thorough  description  of  the  various  noise  sources,  the 
reader  is  referred  to  the  articles  by  Lowenstein  (7)  and  Sakai  (g),  but  the 
effects  of  two  noise  sources  (photon  and  digitization)  in  the  Lincoln  Labora¬ 
tory  experiment  are  of  immediate  interest. 

Where  intense  sources  are  combined  with  large  optical  throughputs,  as  in 
the  case  under  discussion,  the  system  may  become  background  noise  limited 
(photon  or  related  noise  from  the  high  temperature  mercury  arc  lamp).  For 
example,  with  an  incident  total  continuum  power  of  10  ^  W,  the  lamp  output 

g 

stability  must  be  better  than  1  part  in  10  for  one  second  durations,  in  order 
to  insure  that  the  noise  is  less  than  the  Ge  bolometer  detector  noise  of 
5  x  10'13  W. 

Another  serious  noise  contribution  arises  from  the  limited  resolution 
capability  of  FTS  instruments  (see  Fig.  10  showing  a  water  vapor  FTS  spectrum 
taken  over  a  2-m  path  at  17  Torr).  If  the  maximum  resolution  fraction  5v/Av 

(the  ratio  of  the  spectral  linewidth  to  the  resolution  element  width)  is  less 
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than  10  ,  the  absorption  signal  will  appear  m  the  third  figure  of  the  digital 
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output.  This  means  that  the  dynamic  range  of  the  analog-to-digital  conversion 
system  must  be  accurate  to  at  least  four  figures,  or  fluctuations  in  the  third 
digit  will  produce  "digitization  noise"  that  will  contribute  sufficiently  that 
the  multiplexing  advantage  will  be  lost  and  the  sensitivity  reduced  accord¬ 
ingly. 

As  discussed  in  the  Appendix,  incoherent  detection  schemes  are  mainly 
suited  to  survey  purposes,  where  high  resolution  is  not  a  requirement.  For 
studying  radiative  transfer  in  gaseous  media,  FTS  is  normally  used  only 
where  collision  broadening  produces  linewidths  large  enough  to  fall  within  the 
resolution  capabilities  of  the  instrument.  Where  Doppler  broadening  is  dominant, 
unresolved  absorption  signals  can  be  diminished  by  orders  of  magnitude,  and 
the  requirement  for  high  dynamic  range  of  analog-to-digital  converters  add 
further  complications  to  the  detection  problem,  because  the  absorption  data 
may  be  contained  only  in  the  fourth  or  fifth  digit.  A  spectrometer  obtained 
commercially,  such  as  the  one  available  for  the  plume  experiments,  may  require 

expensive  modifications  to  increase  the  dynamic  range.  When  a  single  spectral 
line  is  the  object  of  investigation,  a  high-resolution  heterodyne  receiver  is 

clearly  a  more  suitable  instrument,  particularly  at  lower  infrared  frequencies. 

C.  Heterodyne  Detection 

The  theoretical  sensitivity  of  a  heterodyne  receiver  operating  in  the 
Rayleigh- Jeans  region  (in  one  sense  of  polarization)  may  be  deduced  from 
Eq.  ( A- 10 ) : 


_  kAT  MIN(  6v,B.„) 

(  )  - 

n  het 


(26) 
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where  is  the  antenna  efficiency,  including  the  beam-fill  factor,  and  n,, 
is  again  the  observation  efficiency.  MIN(5v,  B^)  denotes  the  smaller  of  the 
two  quantities:  <5v  and  B^.  The  above  relation  may  be  simplified  further  by 
assuming  that  the  spectral  linewidth  6v  =  B^,  the  IF  (intermediate  frequency) 
bandwidth.  This  condition  will  provide  the  maximum  signal-to-noise  ratio  be¬ 
cause  there  is  no  loss  of  signal  per  resolution  width  from  lack  of 
resolution,  which  occurs  when  increases  beyond  <Sv,  and  no  further  increase 
in  noise  per  width  from  the  existing  /b_^  dependence,  which  occurs  when 
decreases  below  6v.  Although  this  latter  case  would  provide  more  details 
of  the  line  structure,  it  would  do  so  only  at  the  expense  of  reduced  sensitivity 
{ see  Appendix ) . 

In  previous  work  (l),  it  was  shown  that  6v  ~  10~^v  for  Doppler-broadened 
spectral  lines,  and  for  single-polarization,  single- sideband  receivers, 

(NEP),  ,  ~  hv/n  ,  where  n  is  the  quantum  efficiency.  With  these  substitutions 
for  <$\)  and  (NEP)^^,  Eq.  (26)  may  be  written  as 

(|)  =  n,  n  /ru(l-n.,)  2.1  x  io3  AT  /£  ,  (27) 

N  het  1  q  2  v 

where  v  is  given  in  GHz.  For  comparison  with  the  FTS  instrument,  consider  the 
following  parameter  values: 

v  =  752  GHz 

AT  -  T.  -  T  =  3000  K  -  50  K  =  2950  K 
b  p 
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t  =1  sec. 


ru  =0.3  (assuming  100$  beam  fill) 

nq  =  5.0  x  io"3  (i.e.,  (NEP)het  =  10-19  W/Hz) 

n2  =0.5 

Under  these  laboratory  conditions,  Eq.  (27)  yields 

(§)  =  169  (28) 
het 

provided  that  the  receiver  is  detector-noise  limited.  (For  the  satellite 
radiometer  case,  =  250  K  and  (S/N)het  =  14.) 

At  submillimeter  wavelengths,  the  heterodyne  approach  offers  definite 
advantages  over  the  FTS  instrument,  in  the  situation  where  only  a  single  spec¬ 
tral  line  is  to  be  investigated.  In  practice,  the  theoretical  sensitivity  of 
a  coherent  receiver  is  more  easily  realized  because  the  resolution  requirements 
are  met  by  adjusting  EL^  —  5v.  As  a  result,  large  dynamic  ranges  of  readout 
instruments  are  not  required  and  the  problem  of  digitization  noise  is  not  en¬ 
countered. 

However,  to  construct  a  heterodyne  receiver  for  operation  at  submilli¬ 
meter  wavelengths,  there  are  two  difficult  technological  requirements  which 

must  be  met:  1)  a  mixer  with  low  (NEP)^^  values,  in  this  case  about 
-19 

10  W/Hz  for  557-GHz  water  lines,  and  2)  a  local  oscillator  (L0)  capable 
of  delivering  milliwatts  of  power  at  frequencies  close  enough  to  the  spectral 
line  of  interest  for  the  use  of  microwave  IF  frequencies. 
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In  Fig.  11,  a  sketch  of  a  radiometric  apparatus  using  a  receiver  of 
-19 

(NEP)  os*  10  W/Hz  is  shown  for  the  measurement  of  a  77  K  blackbody 
het 

(immersed  in  liquid  nitrogen)  (9).  The  input  signal  is  first  chopped  by  a 
room  temperature  blade  and  then  merged  with  a  local  oscillator  (L0)  signal 
from  a  submillimeter-wave  laser  through  a  diplexer.  By  means  of  an  ellip¬ 
soidal  or  paraboloidal  mirror,  the  input  and  L0  signals  are  focused  onto  a 
Schottky-diode  mixer  mounted  in  a  quasi-optical  antenna  fixture  referred  to 
as  a  comer-cube  reflector  (9).  The  details  of  the  diode  and  reflector 
structure  are  shown  in  Fig.  12  with  its  antenna  pattern  data.  Local  oscillator 
power  can  be  obtained  from  the  submillimeter-wave  laser  sketched  in  Fig.  13, 
where  a  40- W  C02  laser  delivering  a  10. 6- pm  wavelength  beam  is  used  to  furnish 
optical  pumping  of  a  far-infrared  (FIR)  laser,  an  end-mirrored  cell  which  con¬ 
tains  gas  at  a  pressure  on  the  order  of  200  mTorr.  To  observe  absorption  by 
any  water-vapor  rotational  line,  this  gas  must  be  selected  to  provide  FIR  laser 
output  at  a  neighboring  frequency.  For  the  557-GHz  line,  difluoroethane 
2 ^  produce  power  at  563  GHz;  for  the  752-GHz  line,  formic  acid 
(HC00H)  will  deliver  a  strong  output  (30  mW)  at  761  GHz. 

An  attractive  alternative  to  the  submillimeter-wave  laser  L0  is  to  use 
a  harmonic  frequency  of  a  carcinotron  power  tube.  For  example,  a  carcinotron 
producing  1  watt  of  output  power  at  280  GHz  may  be  employed  for  second  har¬ 
monic  (560  GHz)  generation  of  several  milliwatts,  more  than  adequate  to  serve 
as  the  L0  for  a  557-GHz  receiver.  Harmonic  mixing  is  a  powerful  tool  in  milli¬ 
meter  and  submillimeter-wave  technology  and  may  also  be  used  effectively  for 
converting  the  heterodyne  receiver  to  a  spectrometer. 
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In  Fig.  14,  a  second  Schottky  diode  is  employed  to  convert  the  heterodyne 
receiver  into  a  spectrometer  by  first  mixing  the  LO  signal  with  an  inde¬ 
pendently  generated  IF  signal  to  produce  two  sidebands  at  frequencies 

and  -  v^.  In  this  manner,  three  submillimeter  signals  may  be 
radiated  through  the  plume,  one  of  which  is  tuned  to  the  desired  water  line 
frequency  by  adjusting  v.^..  At  the  opposite  port,  the  second  mixer  may  be 
used  to  recover  the  IF  signal  from  the  two  sidebands,  and  produce  a  very 
sensitive  and  highly  selective  ( determined  by  the  narrow  bandwidth  of  the 
initial  IF  signal)  measurement  of  the  plume  absorption.  By  sweeping  the  IF, 
a  continuous  display  of  the  absorption  line  may  be  established  (10  ). 

Although  this  type  of  spectroscopic  measurement  will  not  by  itself 
represent  the  basic  radiometric  situation  described  previously  ( 1 ),  it  is 
a  far  more  powerful  technique  for  plume  laboratory  diagnosis,  because  it  not 
only  produces  a  more  intense  background  signal  ( the  effective  temperature  has 
been  estimated  at  1010  K),  but  also  offers  the  capability  of  spatial  mapping, 
since  the  off-axis  paraboloid  used  to  collimate  the  signal  will  produce  a 
beam  only  a  few  centimeters  in  diameter  at  the  center  of  the  plume. 
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Fig.lL.  Submillimeter-wave  spectrometer  for  high  sensitivity,  high 
frequency  resolution  measurements  of  plume  water  vapor  absorption 
line. 
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IV.  CONCLUSIONS 


Based  on  the  best  theoretical  estimates  calculated  from  the  PLUME  1  com¬ 
puter  program,  cold  rocket  plume  core  signatures  thousands  of  nozzle  diameters 
in  length  should  exist  for  the  557-GHz  and  752-GHz  lines  of  water  vapor.  Sig¬ 
natures  of  this  magnitude  should  make  it  possible  for  third  stage  and  mid-course 
vernier  motors  a+  250  km  altitude,  or  bus  motors  at  350  km,  to  be  observed  with 
antenna  apertures  of  several  meters  at  dwell  times  of  one  second.  Because  of 
the  transparency  of  the  high  temperature  mixing  layer  which  surrounds  the 
cold  core,  the  contribution  from  this  region  is  not  expected  to  have  any  sig¬ 
nificant  impact  on  the  effective  temperature  of  the  signature.  For  the  case 
of  the  core  optical  depth  on  the  order  of  unity,  the  mixing  layer  optical 
depth  is  about  20  times  smaller. 

With  the  recent  development  of  a  submillimeter-wave  mixer  capable  of  pro- 

-19 

ducing  heterodyne  receiver  noise  equivalent  powers  of  10  W/Hz,  a  radiometer 
for  detecting  cold  plumes  against  a  warm  earth  background  is  a  technological 
reality.  For  a  background  temperature  of  250  K,  the  cold  core  should  be  ob¬ 
servable  with  a  signal-to-noise  ratio  of  14  from  synchronous  orbit  with  post¬ 
detection  integration  times  on  the  order  of  1  second.  Dwell  times  of  only  10  ms 

-20 

should  become  possible  when  mixers  with  NEP  values  on  the  order  of  10  W/Hz 

become  available.  Self-emission  from  the  plume  core  is  also  detectable,  but, 

for  the  same  signal-to-noise  ratio,  only  with  dwell  times  that  are  longer  by 

the  square  of  the  ratio  AT/T  ,  where  AT  —  -200  to  -150  K  when  T  ~  50  to  100  K, 

P  p 

depending  on  plume  altitude. 
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For  laboratory-scale  simulations,  the  requirements  for  beam  fill  are  not 
difficult  to  meet  because  of  the  proximity  of  source  to  antenna.  In  addition, 
FIR  lasers  or  carcinotrons  with  harmonic  mixing,  both  cumbersome  and  inefficient 
LO  power  sources,  may  be  used  with  relative  convenience  in  a  laboratory  environ¬ 
ment.  A  laboratory  demonstration  of  the  radiometric  capabilities  of  the  hetero¬ 
dyne  receiver  for  observing  cold  rocket  plumes  at  557  and  752  GHz  should  be 
possible  with  these  recent  technological  developments. 

For  field  measurement,  the  extreme  opacity  of  the  lower  atmosphere 
(troposphere)  at  557  and  752  GHz  renders  unlikely  the  observation  of  rocket 
plumes  from  an  aircraft  at  these  frequencies.  Only  by  means  of  a  space  plat¬ 
form,  e.g.,  the  Space  Shuttle  or  high-altitude  rocket  probe,  may  a  proper  field 
measurement  be  carried  out,  and  only  at  altitudes  high  enough  to  be  capable 
of  viewing  the  target  against  a  250-K  earth  background  or  against  the  cold 
sky.  Therefore,  field  demonstrations  of  the  feasibility  of  the  submillimeter- 

wave  cold  plume  detection  concept  will  be  dependent  on:  1)  the  availability  of 
a  suitable  space  platform,  and  2)  the  eventual  development  of  an  efficient, 
compact  and  preferably  tunable  source  of  local  oscillator  power.  Although 
carcinotrons  may  be  operated  from  a  manned  space  laboratory,  this  latter  re¬ 
quirement  strongly  suggests  the  need  for  improved  solid-state  technology 
in  this  area. 
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APPENDIX 


A.  Infrared  Detector  Concepts 

For  passive  detection  of  a  radiating  source  at  temperature  T,  the 
brightness  spectral  distribution  is  given  by  the  Planck  formula 


b( v,T )  = 


2hv  ,  1 

2  1  hu/kT  n  ' 

c  e  -1 


(A-l ) 


-2  -1  -1 

in  units  of  '.'/-cm  -ster  -Hz  ,  where  h  is  Planck's  constant,  k  is  Boltzmann's 
constant,  c  is  the  speed  of  light,  and  v  is  the  frequency. 


When  a  radiating  plume  source  at  excitation  temperature  Tp  is  detected 
against  an  equivalent  blackbody  continuum  background  of  temperature  T^,  the 
resultant  brightness  is 


bp(v,Tp)  =  b(  v,Tp)(  1-e  V)  (A-2) 

bb( V,Tb  )  =  b^v’Tb)’ 


where  is  the  optical  depth  of  the  plume.  The  second  factor,  in  parentheses, 
accounts  for  thin  sources  radiating  less  than  the  effective  blackbody.  The 
resultant  brightness  contrast  may  be  represented  by 


-T 

Ab(v,Tp,Tb)  =  b  (v,T  )(l-e  V) 


+  Vv-Tb)e  Tv  -  V^V  ’ 


(A- 3 ) 
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which  reduces  to 


Ab(v,Tp,Tb)  Mb  (v,T  )  -  bb(v,Tb)](l-e  V)  .  (A-4) 

For  a  source  of  area  A,  which  emits  over  a  solid  angle  n  that  can  reach 
the  detector,  the  signal  power  over  a  bandwidth  Av  is  given  by 

S  =  AbAvAft,  ( A-5  ) 

assuming  100%  collection  efficiency  to  the  detector. 

The  product  Ail  is  commonly  r  .'erred  to  as  the  "throughput"  or  "etendue" 
of  a’.  ;  rical  system  and  is  an  important  factor  in  determining  the  radiation 

gathering  capability  of  a  receiver.  For  incoherent  detection  systems,  A 12  is 
independent  of  frequency  and  is  limited  only  by  the  geometrical  factors  associa¬ 
ted  with  the  practical  size  of  the  source  and  receiver  aperture  (see  Fig.  15). 
With  coherent  or  heterodyne  systems,  the  received  signal  is  that  contained 
within  the  main  antenna  beam  of  solid  angle  fi,  corresponding  to  an  effective 
antenna  area  A,  and  the  throughput 


v 


according  to  the  conditions  of  diffraction  (11). 
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Fig. 15.  Geometrical  optics  definition  of  "throughput"  AS  1.  Quantities 
with  subscript  m  refer  to  the  collecting  mirror,  whose  effective  aper¬ 
ture  diameter  is  M  and  focal  length  is  F,  while  quantities  with  sub¬ 
script  s  refer  to  the  source  (arc  lamp,  e.g.)  of  effective  diameter  S. 
Note  that  the  throughput  for  the  source  equals  that  for  the  mirror: 


A  n 
s  s 


=  A  ft  = 
m  m 
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In  general,  the  rms  noise  of  an  incoherent  system  may  be  expressed  as 


N. 

me 


(NEP).  //t  +  N' 
me 


( A-7 ) 


-1/2 

where  (NEP)^^^  is  the  noise-equivalent  power  in  W-Hz  of  the  detector,  t 
is  the  postdetection  integration  time,  and  N'  is  the  power  contribution  from 
non-random  noise  sources. 

The  corresponding  relation  for  the  rms  noise  of  a  heterodyne  receiver 
is 


Mhet=<"EP)het/h7^/r*"'>  (A-' 

where  B._  is  the  IF  bandwidth  and  (NEP),  .  is  the  noise  equivalent  power  in 
xi  het 

units  of  V/-Hz-^. 

After  combining  Eqs.  (A-5)  through  (A-8)  and  introducing  efficiency 
factors,  sensitivity  relations  may  be  written  as 

^  _ Ab6v  Afl _  (A-1 

N  inc  1  (NEP)inc/An2(i-n2)  +  N' 

and 


,Ss  _  Ab6v  ( c2/v2 ) 

'  M  /  n«i  _  . - — -  -  -  - ~ 

het  1  (NEP)  ^Bjf/tn2(l-n2)  +  N 1 


(A-10) 
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A 


when  exceeds  6v.  When  B.^,  <  6v,  the  spectral  linewidth,  the  total  signal 
in  the  channel  of  width  B^  is  reduced  by  the  fraction  B^/6v  and  represents 
the  condition  of  resolving  the  spectral  line.  When  B^  >  6v,  the  line  is  un¬ 
resolved  and  the  noise  in  the  channel  of  width  B^  is  larger  than  it  need  be. 
Obviously,  the  optimum  signal-to-noise  ratio  occurs  when  B..^~  6v.  The 
efficiency  is  for  the  antenna,  including  the  beamfill  factor,  and  is  the 
fraction  of  the  observation  or  integration  time  t  that  is  devoted  to  the  plume, 
as  opposed  to  the  comparison  or  reference  source  during  the  remaining  time. 

B.  Coherent  Versus  Incoherent  Detection 

For  infrared  detection,  the  relative  merits  of  coherent  and  incoherent 
systems  depend  to  a  large  extent  on  the  frequency  range  and  degree  of  frequency 
resolution  involved  in  the  particular  application.  Basically,  there  are  three 
areas  where  these  schemes  differ,  1)  the  collected  signal,  2)  the  detector  noise, 
and  3)  the  resolution  capability. 

As  described  earlier,  the  amount  of  detectable  signal  power  is  determined 
by  the  throughput  A ft  in  either  case.  For  the  incoherent  system,  A SI  is  indepen¬ 
dent  of  wavelength  and  is  limited  only  by  the  geometrical  factors  associated 

with  the  collection  optics  or  antennae.  With  diffraction-limited  heterodyne 
2 

receivers,  Afi  =  A  ,  as  determined  by  the  main  beam  of  the  antenna,  with  the 
result  that  the  signal  received  is  usually  much  less  than  that  of  an  incoherent 
system,  when  considering  shorter  wavelengths.  Only  in  the  millimeter-wave  or 
far-infrared  regions  will  the  coherent  detector  become  competitive  if  spectral 
resolution  is  not  a  consideration.  However,  because  the  ultimate  NEP  of  a 
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heterodyne  receiver  is  approximately  the  energy  of  a  single  photon  hv,  i.e.,  the 
"quantum  limit"  (12  ),  the  minimum  detectable  signal  power  is  usually  much  smaller 
than  that  of  an  incoherent  system,  particularly  where  the  predetection  band¬ 
width  B.^,  is  small.  This  feature  of  heterodyne  detection  offsets  its  lower 
signal -collecting  capability  to  some  extent  and  further  enhances  its  sensi¬ 
tivity  at  lower  frequencies. 

Finally,  for  high  frequency-resolution  applications,  heterodyne  receivers 
have  a  considerable  advantage  over  incoherent  systems.  For  diffraction  grating 
and  Fourier  transform  spectrometers,  the  minimum  resolution  cell  Av  is  on  the 
order  of  a  few  GHz,  while  the  resolution  limit  of  a  heterodyne  receiver  is  de¬ 
termined  by  the  values  of  B^,  which  may  be  less  than  a  MHz.  Because  Av  for 
a  grating  instrument  is  determined  by  the  aperture  slit  width,  the  magnitude 
of  A  is  proportional  to  Av  and  it  may  be  argued  from  Eq.  (A-9)  with  Av  =  6v, 
the  linewidth, 


(f)  -  («v)2  ,  ( A-ll ) 

grat 

while  for  an  FTS  the  aperture  is  fixed,  and 


(S) 


6v 


FTS 


(A-12) 


However,  since  =  6v  for  maximum  sensitivity  of  a  heterodyne  receiver, 
Eq.  (A-10)  indicates  that 
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(-) 


het 


(A-13) 


Thus,  the  heterodyne  detector  is  not  only  capable  of  higher  resolution,  but 
its  sensitivity  is  less  dependent  on  6v  and  will  not  decrease  as  severely 
when  high  resolution  is  achieved. 
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